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Abstract We review here some magnetic phenomena in astrophysical parti- 
cle accelerators associated with coUisionless shocks in supernova remnants, 
radio galaxies and clusters of galaxies. A specific feature is that the accel- 
erated particles can play an important role in magnetic field evolution in 
the objects. In particular, we discuss a number of cosmic-ray (CR) driven, 
magnetic field amplification processes that are likely to operate when diffu- 
sive shock acceleration (DSA) becomes efficient and nonlinear. The turbulent 
magnetic fields produced by these processes determine the maximum ener- 
gies of accelerated particles and result in specific features in the observed 
photon radiation of the sources. Equally important, magnetic field amplifi- 
cation by the CR currents and pressure anisotropics may affect the shocked 
gas temperatures and compression, both in the shock precursor and in the 
downstream flow, if the shock is an efficient CR accelerator. Strong fluctua- 
tions of the magnetic fleld on scales above the radiation formation length in 
the shock vicinity result in intermittent structures observable in synchrotron 
emission images. The flnitc size twinkling, intermittent structures - dots, 
clumps, and filaments - are most apparent in the cut-off region of the syn- 
chrotron spectrum. Even though these X-ray synchrotron structures result 
from turbulent magnetic fields, they could still be highly polarized provid- 
ing an important diagnostic of the spectrum of the turbulence. We discuss 
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both the thermal and non-thermal observational consequences of magnetic 
field amplification in supernova remnants and radio-galaxies. Resonant and 
non-resonant CR streaming instabilities in the shock precursor can generate 
mesoscale magnetic fields with scale-sizes comparable to supernova remnants 
and even superbubbles. This opens the possibility that magnetic fields in the 
earliest galaxies were produced by the first generation Population III super- 
nova remnants and by clustered supernovae in star forming regions. 

Keywords radiation mechanisms: non-thermal — X-rays: ISM — (ISM:) 
supernova remnants — clusters of galaxies-shock waves-magnetic fields 



1 Introduction 

Particle acceleration takes place in many active astrophysical objects of very 
different nature and scales. Magnetic fields play the central role in charged 
particle acceleration either as a intermediary between the plasma flows and 
energetic particles, as is the case for Fermi acceleration (e.g., in coUisionless 
shocks), or as a source of free energy to be converted into energetic particles 
(i.e., magnetic field reconnection processes). 

The existence of the highly amplified magnetic fields in the range of 0.1-1 
mG in the shells of young supernova remnants (SNRs) was established as- 
suming equipartition between relativistic pa rticles and magnetic fields m the 
synchrotron radio emitting shells (see e.g.. iGinzburg and Svrovatskiil [l963. 
and the references therein). The topology of the magnetic field in SNRs in- 
ferred from the observations of synchrotron emission differ strongly between 
young and old SNRs. Radio polarization studies reveal super-adiabatic mag- 
netic field amplification and a net radial orientation of the magnetic fields in 
young SNRs, while it i s often just shocked interstellar field, mainly tangen- 
tial, in old SNRs (e.g.. lMilnjll990l) . Milne has also pointed out that as the 
resolution increases, the polarization structure becomes more complex. 

A pixel- by- pixel map of Faraday rotation has been produced applying 
rotation measure (RM) synthesis to the data observed with the Australia 
Te lescope Compact Ar r ay, fo r the entire supernova remnant G296. 5-1-10.0 
bv lHarvev-Smith et al.l (l2010l ). A highly ordered rotation measure structure, 
with an anti-symmetric rotation measure pattern, was observed. The authors 
propose that the observed rotation measures are the imprint of an azimuthal 
magnetic field in the stellar wind of the progenitor star. A swept-up magne- 
tized wind from a red supergiant can produce an azimuthal pattern of the 
magnetic field at large distances from the star and can naturally produce the 
observed anti-symmetric RM pattern. Supernova expansion into such a wind 
could account for the apparent bilateral structure of the SNR's radio and X- 
ray morphologies. In the case of SN1006, a comparison between observed and 
synthesized radio maps, making different assumptions about th e dependence 
of elec tron injection efficiency on the shock obliquity, allowed iPetruk et al.l 
( 20091 ) to constrain possible nonthermal electron injection models. 

A number of possible amplification mechanisms were considered. Rayleigh- 
Taylor instabilities in the shell have been considered for many years as a 
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potential source of the turbulent magnetic fields. The Rayleigh- Taylor insta- 
bility at the interface of the ejecta an d the shocked ambient m edium was 
proposed to explain these observations. lJun and Normaiil (|l996f ) performed 
multi-dimensional MHD simulations of the instability in the shell and its 
effect on the local magnetic field. They found that the evolution of the in- 
stability is very sensitive to the deceleration of the ejecta and the evolu- 
tionary stage of the remnant. The Rayleigh- Taylor and Kelvin-Helmholtz 
instabilities amplify ambient magnetic fields in the simulations locally by as 
much as a factor of 60 around dense fingers due to stretching, winding, and 
compression. Globally, the amount of magnetic-field amplification was nev- 
ertheless low and the magnetic energy density reaches only about 0.3% of 
the turbulent energy density at the end of simulation. Strong magnetic field 
lines draped around the fingers produce the radial B-vcctor polarization. 
M agnetic Rayleigh- Tay l or ins tabilities in three dimensions were simulated 
by IStone and Gardiner! (|2007l ). with a focus on the nonlinear structure and 



evolution that results from different initial field configurations. They found 
that strong magnetic fields do not suppress the instability but reduce mixing 
between the heavy and light fluids and cause the rate of growth of bubbles 
and fingers to increase in comparison to hydrodynamics. ISchure et all (|2009[ ) 
suggested that higher compression ratios and additional turbulence due to 
the dominant presence of CRs can be important in attempts to reproduce 
the observed magnetic field morphology with the Rayleigh- Taylor instability 
simulations. Both numerical simulations and laboratory laser experiments 
are ongoing to study the RayleighTaylor instability problem. The possible 
effects of magnetic fields on laser experiments of Ray leighTaylor instabil i- 
ties in supernova-like setups was discussed recently bv iFrvxell et al.l ( 2010l) . 



Comparing numerical simulations, as well as experimentally using the Omega 
Laser at the University of Rochester, the authors discussed a possibility that 
the morphology of the instability is significantly altered by the generation of 
very strong magnetic fields during the laser experiments. 



On the other hand. IChevalieil ( 1977t ) shown that radio observations of Ty- 



cho's supernova remnant indicated the presence of a coUisionless shock wave 
undergoing turbulent magnetic field amplification by a factor of about 20. 
He noted that the amplification resembled some phenomena in heliospheric 
coUisionless shocks. In fact, direct support for the sel f-generation of MHD 
turbu lence has been seen in heliospheric shocks, e.g., iBamert et al ] (I2004L 



I2OO8I ). In the upstream region of a shock associated with the Bastille Day 
coronal mass ejection, the excitation of hydromagnetic waves with power 
spectral den sity levels well above the levels in the ambient solar wind were 
observed bv lBamert et al ] (I2004D_ with the magnetome ter on board the Ad- 



vanced Composition Explorer. Achterberg et afl (|1994[ ) set upper limits on 



the scattering mean free paths of radio emitting electrons in front of super- 
nova remnant shock fronts using high-resolution radio observations of four 
Galactic supernova remnants. They found that, for the sharpest synchrotron 
radio rims, the mean free path is typically less than one percent of that de- 
rived for cosmic rays of similar rigidity in the interstellar medium, implying 
that enhanced hydromagnetic wave intensity is generated by diffusive shock 
acceleration. 
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Recent observations of X-ray synchrotron radiation from several young 
SNRs revealed thin X-ray filaments in the vicinity of blast waves. This pro- 
vides indirect, but rather strong evidence for strong super-adiabatic magnetic 
field amp lification (MFA) in the blast wave associated with CR production 
(see e.g., Vink and Laming '2003', 'Bamba et al.' '2005', '2006|, lUchivama et al.l 
l2007LlPatnaude and Fesen2008 , Rcvnolds 2 008 , Vjnk 200a) . The synchrotron 
filamentary structures could be due to a narrow spatial extension of the high- 
est energy electron population in the shock acceleration region limited by 
rapid synchrotron energy losses. In this case, strong magnetic field amplifica- 
tion in the shock vicinity is required to match the observed thickness of the 
non-thermal filaments (e.g.. lVink and Laminej[2b03llBamba et al.ll2005( ). thus 
supporting the case for highly efficient DSA. An alternative interpretation, 
that the observed narrow filaments are limited by magnetic field damping and 
not by the energy losses of the radiati ng electrons, h as also been proposed 
(see iPohl et al.l 120051 l&att et al.ll2008l for example). [C assam-Chenai et al.l 
(|2007l . l2008[) lmve studied this effect in detail for Tycho's SNR and SN1006. 

Diffusive shock acceleration is a very promising mechanism for produc- 
ing superthermal and relativistic particles in a wide r ange of astrophysi 
cal objects ranging from th e Earth bow shock (e.g., Ellison et al.l Il99 



to shocks in galaxy clusters (iBlandford and Eichleilll987 . Jones and EUiso 



ll99lL iMalkov and Drurvll200ll iBvkov et al.ll2008a^ . This mechanism is be- 



lieved to be efficient (e.g.. iHelder et al.l2009| ) and capa ble of producing CR s 
to energies well above 10"'^^ eV in young SNRs (e.g.. iPtuskin et al.ll2010D . 
and even higher in active radio-galaxies such as Cen A (see, for example, 
ICroston et al]|2009[ ). 

Fast and efficient DSA demands that particles are multiply scattered by 
magnetic fluctuations in the shock vicinity. The amplitude of the required 
MHD turbulence is substantially higher than the ambient MHD turbulence 
suggesting that the turbulence is generated by the shock. In order for DSA 
to accelerate particles to high energies, the energetic particles must be able 
to interact with magnetic turbulence over a broad wavelength range. The 
weakly anisotropic distribution of accelerated particles, is believed capable 
of producing this turbulence in a symbiotic relationship where the magnetic 
turbulence required to accelerate the CRs is created by the accelerated CRs 
themselves. In efficient DSA, this wave-particle interaction can be strongly 
nonlinear where CRs modify the plasma flow and the speciflc mechanisms of 
magnetic field amplification. 



2 Magnetic field amplification by CR-driven instabilities 

The non-equilibrium distributions of accelerated particles in efficient DSA 
arc believed capable of producing magnetic turbulence and magnetic field 
amplification. For the forward shock in a young, isolated SNR, the acceler- 
ated CR pressure can be a sizeable fraction of the shock ram pressure and far 
exceed the background plasma pressure. The anisotropic distribution of CRs 
in the shock precursor results in ponderomotive forces on the background 
plasma that can cause the amplification of magnetic fluctuations in a certain 
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wavevector range and polarization. We briefly discuss in section 12.11 the ki- 
netic background for the description of this CR anisotropy and then outline 
a number of instabilities related to the anisotropic CR distributions. 



2.1 Anisotropic CR distribution in DSA 

First-order Fermi acceleration (also called diffusive shock acceleration) of CRs 
by non-relativistic magnetized flows, is characterized by a nearly isotropic CR 
momentum distribution in a broad energy range. Two notable exceptions to 
this are injected, barely superthermal particles, and the highest energies CRs. 
The CR distribution is formed by the multiple scattering of charged parti- 
cles across the shock by magnetic fluctuations over a wide dynamic range. 
Cosmic-ray particle dynamics in DSA is governed by electro-magnetic fields 
both regular and stochastic. The distribution function of energetic particles, 
/(r,p,t) = {F{r,p,t))c, averaged over electro-magnetic fluctuations with 
scales below some value Ic, appropriate for a particular problem, satisfies the 
kinetic equation 

where O is the momentum rotation operator, defined by 

£ is the particle energy, and E = (Eo)c and B = (Bo)c are the electric and 
magnetic fields averaged over electro-magnetic fluctuations with scales below 
Ic- The fluctuating parts of the electromagnetic fields are b = Bq — (Bo)c and 
e = Eg — (Eo)c- The fields are self-consistent since the sources in Maxwell's 
equations include the CR currents determined by F{r, p, t). We should stress 
here that it is not trivial to get a closed system of equations for the averaged 
distribution function /(r, p, t) = (F(r, p, t))c since the higher order moments 
of F may contribute through {h{r,t)F{r,p,t))c and {e{r,t)F{r,p,t))c. In 
fact, even the existence of a well defined averaging procedure, and the aver- 
ages themselves, is assumed rather than proved. 
The collision operator 

I [/] = -Ze{e ~ - |b . dF),, (3) 

is, in general, a nonlinear operator if the effect of energetic particles on the 
electromagnetic fields is non- negligible. 

Different approaches are used to get a closed analytic expression for Ic [f] 
including Fokker-Planck type equations and quasi-linear theory (see e.g., 
Toptvginlll985llBlandford and Eichlejl987llBCTezinskn et al.ll99(]llSchhckeiseij 



The simplest case for the collision operator can be deduced from 



Eq. for particles scattered by magnetic fluctuations with scales, Ic, smaller 
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than the particle gyroradius, rg = cp/ eB, in the total magnetic field B (ran- 
dom plus coherent). In the high-energy limit, Tg » Ic, the particle momen- 
tum change over the field correlation length Ic is small (oc Ic/tq) and the 
scattering rate in Eq. ([3|) reduces to v/2A{p) with A{p) — gRgj(p)/lc. 
Here, Rstip) is the particle gyroradius in the stochastic magnetic field, and g 
is a numerical factor depending on the short-scale fluctuation spectral index. 
The strong momentum dependence, A cc p^, makes this scenario unfavorable 
for confining and accelerating high-energy particles in DSA. 

It is important to note that can be reduced to the angular part of 
the Laplace operator with well-known spherical harmonics as the eigenfunc- 
tions. The high-energy asymptotic range is valid even in the case of strong 
(but quasi-static) magnetic fluctuations of amplitudes larger than the mean 
magnetic fleld. However, for particles of intermediate and low energies, with 
Tg < Ic and scattered by strong magnetic fluctuations, the collision oper- 
ator is no longer simple and only some phenomenological scheme like the 
relaxation time approximation can be used. These schemes do not have a 
microscopic justification. 

The momentum distribution of accelerated particles, /(r, p, t), of a weakly 
anisotropic distribution can be presented as 

7? 

/(p) = ^iV(p)[l + J,(p)], (4) 

where ricr is the CR number density, N{p) is the isotropic and (5a (p) is the 
anisotropic part of the distribution function defined to satisfy 

Js^{p)dQp^O, j f{p)p^dpdnp = n„ . 

Due to CR particle scattering by magnetic fiuctuations of different scales, the 
distribution function is nearly isotropic on scales larger than the CR mean 
free path, A(j)), and, therefore, its anisotropic part is small, i.e., 5a (p) <C 1 on 
scales kAij)) < 1. The actual angular dependence of (5a (p), and the momen- 
tum dependence of A(jj)^ are determined by the exact form of the collision 
operator, the structure of regular fields, and the boundary conditions. It can 
be presented in the spherical harmonic expansion 

'5a(p) = J2 '^l^r.M Yl^rni^, if), \m\ < I, (5) 

l,m 

where •& and ip are the pitch and azimuthal angles of the particle momen- 
tum, correspondingly. Since the isotropic part of the distribution function 
is separated from the anisotropic part in Eq. (|4]), Aq^q = 0. The current- 
type anisotropy that is determined by Z\i,m < 1, is usually considered to 
be the dominant term. That results in the standard diffusion approxima- 
tion and the transport equation for the isotropic part N(r,p,t). In the case 
of a high-energy particle scattering off of small-scale fiuctuations, rg ^ Ic, 
the spherical functions are the eigenfunctions of the collision operator since 
02 Yi „^ = -l{l + l)y,,,„. In that case, indeed, Z\2,m ~ ^im- However, the 
high multipole components of the distribution can be substantial in the case 
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of anisotropic scattering, or in plasma flows with a complex structure of 
quasi- regular magnetic fields (e.g., in the presence of regular or stochastic 
magnetic traps). In the latter case one can expect Z\2,m ^i,m and the 
transport equation for the isotropic part of the distribution function A'^(r, p, t) 
can differ from the standard advection-diffusion equation that is used in the 
analytic DSA models. 

Standard analytic DSA models use a simplified description of particle dif- 
fusion. Generally, the diffusion approximation is made where it is assumed 
that a small anisotropy of the momentum distribution can be completely 
described by the CR current and higher mo ments of the momentum an- 
gular distribution can be neglected (see e.g., jBlandford and Eichled [l987l 



iJones andElhi^ll99lL[Malkov an^ Drurvll200lh . Then, further assuming 
that magnetic fields are frozen-in the plasma, the kinetic Eq. ([1]) can be 
reduced to a transport equation for the isotropic part of the distribution 
function N{r,p,t). This approach can account for particle advection, diffu- 
sion and acceleration by the first-order Fermi mechanism in non-relativistic 
shocks. The diffusion coefficient in these analytical models is parameterized 
as D{p) oc Tg, where the numerical index is generally in the range < a < 2. 
The widely used phenomcnological Bohm diffusion model assumes a = 1, i.e., 
the mean free path A{p) ~ rjrgij)), where the factor i] > 1. 

The CR current and anisotropy in the rest frame of a non-relativistic 
shock is 1^1, ml ^ I's/c « 1, where Vs is the shock velocity. The CR 
anisotropy can also be parameterized by the CR drift velocity Udrift relative 
to the background plasma in the rest frame of the upstream plasma. Various 
parameterizations of the CR current in DSA were used in recent studies of 
magnetic field amplification mechanisms by (e.g.. Belli |2004 iPelletier et al.l 



l2006l . lAmato and Blasill2009l iLuo and Melrosell2009l )" based on the CR trans- 
port equation. Monte Carlo simulations of DSA provide an alternative ap- 
proach to the kinetic theory. In Monte Carlo techniques, particle scattering 
rates are parameterized as a function of energy and then the CR current 
(and other CR anisotropy moments) are derived in nonlinear shock models 
accounting for the CR backreactio n on the plasma flow (e.g., lEUison et al.l 
Il996l . IVladimirov et al ] l2006ll2008ll . 



2.2 Resonant CR-streaming instability 

Resonant interactions have long been k nown to amplify magnetic fluctuations 
on th e scale of the CR gyroradius fe.g.. IWentzelll974 [Z^ii3elll2003l . iKulsrudI 
I2OO5I ) . The simplest case of interest for DSA is Alfven wave amplification by 
streaming CRs. If the streaming speed of the cosmic rays through the ambient 
medium exceeds the Alfven speed, then the amplitude of resonant Alfven 
waves grows exponentially with time. The waves are emitted by the CR 
particles gyrating with gyrofrequency f2 as coherent cyclotron radiation. The 
resonance condition is k v cos(??) = f2. This resonance condition is typical for 
DSA. The growth rate is positive if, in the wave frame, the particles stream in 
the opposite direction to the motion of the background fiuid. Equivalently, the 
waves grow if, in the fluid frame, the mean streaming velocity of the particles 
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exceeds the wave velo c ity, i. e., Wdrift > Vph- The growth rate estimated by 
iBlandford and Eichled (|l987t) is 
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Wdrift 



1 



Vph 



(6) 



where ricr is the total number density of resonant cosmic rays and ria is the 
background ion density. 



2.3 Non-resonant streaming instabilities: short-scale 

In addition to resonant amplification, the CR current can also efficiently 
amplify non-resonant magnetic fluctuations with scales smaller than the CR 
gyroradius ( Belli I20d3 ) . The following two points arc important to produce 



strong, short-scale turbulence. 

(i) The cosmic-ray current has only a weak response to the fluctuations 
with scales shorter than the gyroradii, r^o, of CR particles dominating the CR 
current (short-scale fluctuations), i.e., fcr^o >> 1 (see discussion in section 
12. ip . In this case, the electric current of accelerated particles J''"' is an external 
current in the momentum equation of the background plasma. The current is 
governed by sources of energetic particles and by local electromagnetic fields. 
The current J'^'' initiates a compensating return current in the background 
plasma. 

(ii) The ponderomotive force of the CR current, J'^'^ x B/c, is large enough 
to dominate the magnetic field tension in the momentum equation of the 
background plasma. Then the J'^' x B/c force from the CR current drives the 
Bell short-scale instability at scales below the CR gyroradius. This condition 
is satisfied for scales k < ki where 

C £io 

Note that the condition for mode growth, Vg^k > 1, together with Eq. ([7|), 
implies that the CR current that determines fci must be large enough to 
drive the non-resonant instability. The anisotropy of the relativistic particle 
distribution, |/\i,o|, should exceed the ratio of the mean magnetic field energy 
density to the energetic particle energy density Ecr- That is 

|^''"| > si' 

wh ere t h e CR curren t is given by J" ~ \Aifi\encrC. 

iBelll ( 20041 . 120051 ) discovered that the system is unstable against linear 



perturbations that are oc exp(7i -I- zkr) and the return current drives nearly 
purely growing modes. Here, 7 is the linear growth rate. 

In a cold plasma with sound s peed an, w ell below the Alfven velocity Va, 



i peea an, w e 

the linear growth rate obtained bv lBelll ( 20051 ) depends only on the wavevector 



projection kz on the local mean magnetic field, i.e., 

7 ~ 7ma:^kz/k, (9) 
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where 

7max = VaVki\k\ ~ . (10) 

According to the hncar analysis of lBelll ( 2005[ ). the wavenumber of a grow- 
ing mode must satisfy the condition r~Q < k < ki. Therefore, the instabihty 
growth rate is higher than the Alfven frequency Vak. The amphtudes of the 
growing mode in the hnear approximation can be expressed as 

provided that the kinetic energy density in the growing mode dominates over 
the magnetic energy density because fci > k^. This is in contrast to Alfven 
modes where the energy densities are equal. 

The cold plasma approximation used in the analysis of Bell turbulence is 
valid if 

— V>^i^.o-, (12) 

where VTi is the thermal ion velocity. The approximation is typically good for 
galactic SNRs in a warm interstellar medium where the plasma parameter 
j3 = a^lv i ~ 1.0. In the case of PSA by la rge-scale shocks in superbubbles 
see e.g., iBvkov and ToptvginI [l993l l20Qlh . or in hot intercluster plasma, 
the thermal corrections to the wave dispersion relations can be essential. A 
thorough dis cussion of the effects of a h ot plasma on the short-scale modes 

was given by IZweibel and Evere tt' (2010^. 

P article-in- cell simulations by Ria uelme and Spitkovskvl ( 2009l . l2010[ ) showed 



that the back-reaction of the amplified field on CRs would limit the ampli- 
fication factor of the Bell instability to less than a factor of about 10 in the 
upstream fiows of galactic SNRs. The authors studied the possibility of fur- 
ther amplification driven near shocks by magnetized CRs, whose Larmor radii 
are smaller than the length scale of the field that was previously amplified 
by the Bell instability in the upstream flow. They found that additional am- 
plification can occur due to the CR current perpendicular to the field. The 
maximum amplification of the instability is determined by the disruption 
of the CR current, which happens when CR Larmor radii in the amplified 
field become comparable to the length scale of the instability. Amplification 
factors up to ~ 45 were estimated in that case. 

An important f eature of the instability dr i ven b y the perpendicular cur- 
rent established bv lRiauelme and Spitkovskvl ( 2010[) is the characteristic de- 



pendence of the amplified field on the shock velocity, oc v'^, which con- 
trasts with the Bell instability acting alone where oc . Different scalings 
ca n, in principle, be constraine d in statistical studies of radio SNRs, as done 
bv lBandiera and PetrukI ( 2010l ). but the results are sensitive to other param- 



eters such as ambient density. It is also possible to constrain these scalings 
by comparing radio and X-ray synchrotron lightcurves of individual young 
supernova shells. In an amplified magnetic field, the synchrotron lightcurve 
will differ from the radio if strong magnetic field damping behind the shock 
occurs. 
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2.4 Non-resonant streaming instabilities: long-wavelength 



The nonresonant, short-scale instability introduced bv iBelll ( 2004[ ). and dis- 



cussed in section [^751 is fast and can strongly amplify short-scale magnetic 
field fluctuations. These strong, short-scale fluctuations are important for the 
dynamics of CR modified shocks and their emission properties, however, to 
study the maximum energies of CR particles in DSA one needs to study the 
amplification mechanisms for long- wavelength fiuctuations of scales rgok < 1. 
Specifically, what are the consequences of long- wavelength fluctuations prop- 
agating in highly turbulent plasma with much shorter scale fluctuations? 
Also, contrary to the short-scale case, the response of the magnetized CR 
current on magnetic field fiuctuations is not small in the limit rppfc < 1. 

Using a multi-scale, quasi- linear analvsis. lBvkov et al.l (|201l[ ) showed that 
the presence of turbulence with scales shorter than the CR gyroradius en- 
hances the growth of modes with scales longer than the gyroradius, at least 
for particular polarizations. They used a mean-field approach to average the 
equation of motion and the induction equation over the ensemble of mag- 
netic field oscillations acco unting for the anisotropy of relativistic particles 
on the background plasma. iBvkov et al.l ( 201lh derived the response of the 



magnetized CR current on magnetic field fiuctuations and showed that, in 
the presence of short-scale Bell-type turbulence, long wavelength modes are 
amplified. In general, the growth rates depend on the mode propagation an- 
gle. For modes propagating parallel to the initial magnetic field, the growth 
rate is 

and these modes have the fastest growth rates for the Bohm diffusion regime 
with ?7 ~ 1. 

The propagation angle ^max of the mode of maximum growth for 77 > 1 

is 

cos6'„iax = l/?7, (14) 
and the maximum growth rate at j] ^ 1 is determined by 

(15) 

The above results were obtained assuming Vpijvg ^ 1. In the case un- 
der consideration 7(fc) ^ i^(fc), and therefore from Eq. ([T5]) one may get the 
validity condition in the form \/ ttNb /4\/ fco /k 1, where Ma is the 

Alfvcn Mach number of the shock. The angular dependence of the growth 
rate in the hydrodynamical regim e depends on the d imensionless collision 
strength 77, as is shown in Fig. [Tlof iBvkov et al.l ( 201lh . The growth rates of 



resonant, non-resonant Bell's, and n on-resonant l o ng wa velength instabilities 
are illustrated in Fig. [2] taken from iBvkov et all ( 201lh where the polariza- 



tion, helicity, and angular dependence of the growth rates are calculated for 
obliquely propagating modes for wavelengths both below and above the CR 
mean free path. The long-wavelength growth rates estimated for typical su- 
pernova remnant parameters are sufficiently fast to suggest a fundamental 
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increase in the maximum CR energy a given shock can produce. It should 
be noted that the phase velocity of the long-wavelength modes is a growing 
function of the wavelength and it is larger than the Alfven speed. That may 
steepen the particle spectra at the highest energy end by reducing the effec- 
tive shock compression at high energies. Also, the high phase velocity of the 
modes could increase the role of Fermi-II acceleration. 

An important point is that the short-scale turbulence can influence the 
large-scale dynamics through the ponderomotive forces imposed on the plasma 
by the turbulence and the CR current. To derive the mean ponderomotive 
force one must average the momentum equation over the ensemble of short- 
scale fluctuations. When this is done, it is seen that the ponderomotive force 
due to the CR current response may result in a long-wavelength instability 
in a way somewhat similar to Bell's instability. 

2.5 Non-resonant firehose instability 

The two CR streaming instabilities, resonant and non-resonant, despite be- 
ing rather different physical mechanisms, both rely on the non-zero, first 
moment of the CR distribution, Ai rn-, i-C, the CR current. The anisotropic 
CR distributions, however, have non- vanishing higher order harmonics, Ai^m- 
The well-known example of the instability due to the high-order anisotropy 
is the weakly magnetized thermal background plasma where the gas pressure 
anisotropy can be large. The systems we consider are known to be unstable if 
the anisotropy is large enough providing P|| > P± + B'^/4:Tr (firehose instabil- 
ity) o r Pj_ > P||-|-i?^/47r (mirror instability) fe.g.. lTreumann and BaumioharmI 
I1997I . and referenc es therein), and can be impo r tant in astrophysical shock s 
(see, for example, iBlandford and Eichle"rl I1987L iBlandford and FunkI |2007^ ■ 
Since the growth rate is oc fc[(P|| - Pj_- SV47I•)/p]^/^ the instability oper- 
ates in hot, high /3 plasmas. 

In a strong supernova shock, the far upstream plasma is cold unless the 
SNR is located in a hot superbubble. Then, the free energy of the back- 
ground plasma pressure anisotropy is not an efficient driver of magnetic field 
amplification. However, the CR pressure in the shock precursor close to the 
shock is expected to reach a sizeable fraction of the shock ram pressure, i.e., 
pCR _ Cpavl/2. In this case, the anisotropy of the CR pressure defined by 
A2^m should play a role. In the hydrodynamieal regime, where the wavelength 
27r/A: is longer than the CR mean free path, A, the growth rate of the firehose 
instability can reach 

^f\k) « l^2,o|/2 vs k ^ VC72 ^ . (16) 

In this estimate, we assumed that Z\2,o oc A\ q. However, as we have pointed 
out earlier, in some plasma flows with a complex structure of quasi-regular 
magnetic fields, or in the flows with regular or stochastic magnetic traps, one 
can expect much larger pressure anisotropy with Z\2.o Zii.o and, therefore, 
higher growth rates. The mirror instability would occur with a growth rate 
similar to Eq. (jl6p if the perpendicular component of CR pressure dominates 
over the parallel pressure. 
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2.6 CR-pressure gradient driven acoustic instability 



Drury and collaborators ( Dorfi and Drurvll985llDrurv and Fallelll986[ ) showed 



that acoustic perturbations in a two-fluid plasma with a CR pressure gra- 
dient can grow if the scale of the CR gradient is shorter than D/gq, where 
ao i s the therma l gas sound speed and D is the diffusion coefBcient (see 



also lChalovlll988l IZaiik et al.lll990l iMalkov and prurvll200ll iMalkov et al] 
l2ni(C " It is argued bv IMalkov and Diamondl (l2009t ) that this acoustic insta- 



bility can be important for DSA. The growth rate of the acoustic instability 
can be estimated as 

^ac ^ \X 1 (17) 

2pa ao 

for k > (see, e.g.. IMalkov et al]l201Cll for a recent discussion). Since 

IVP*^^! ~ (vl/iAcao), then fcmin ~ ({W^ Vs/c)/A. Thus, the long wave- 
length fluctuations of kA < 1 that are the most important to increase the 
maximum energies of CRs accelerated by DSA, can grow only if C{M'^ Vg/c) < 
1. This condition is rather restrictive for young SNRs in the warm interstel- 
lar medium. On the other hand, the short-scale fluctuations amplified by the 
instability are important for heating the plasma upstream of the shock and 
ca n contribute to the pond eromotive force discussed in section \2A[ 

iBeresnvak et al.l (200^ proposed a DSA model in which stochastic mag- 



netic fields in the shock precursor are amplified through small-scale dynamo 
effects. The solenoidal velocity perturbations that are required for the dy- 
namo to work are produced in the model through the interactions of the 
pressure gradient of the CR precursor and the density perturbations in the 
inflowing fluid. 



2.7 Instabilities in partially- ionized plasmas 

The Hq, emission from SN1006 and some other SNRs implies the presence 
of a sizable fraction of partially-ionized plasma in the circumstellar medium. 
Using a kinetic desc ription of CRs, comb i ned w ith a fluid description of the 
background plasma, fevkov and Toptvglnl (|2005| ) and lReville et al.l ( 2007t) in- 



vestigated the growth of hydro-magnetic waves driven by CR streaming in 
the partially-ionized precursor environment of a supernova remnant shock. 
In this case, modifications of Ohm's law in the magnetized plasma can be 
important for the instability, as can ion-neutral damping. If the cosmic ray 
acceleration is efficient, large neutral fractions are required to compensate 
the growth of the non-resonant mode. For modest acceleration rates, the ion- 
neutral damping can dominate over the instability of weakly driven modes 
even at modest ionization fractions. In the case of a supernova shock inter- 
acting with a molecular cloud, such as in IC 443, W28, W44, or W51C, the 
low ionization in the ambient medium could limit the maximum energy of 
accelerated particles s ince this is largely determined by the upstream mag- 
netic turbulence (e.g.. lBvkov et al.ll2000t) . Observationally, SNRs interacting 
with molecular clouds arc indeed found to be sources of high-energy gamma- 
ray emission (see e.g.. (Albert et al..,2007. .Aharonian et al.,,2008^ ,Abdo et all 
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Fig. 1 The angular dependence of the growth rates of the long-wavelength unstable 
modes in the hydrodyamical regime where kvgo < rj~^ , for rj = 10. The angle 6 is 
between the wavevector and the CR current in the upstream rest frame of a parallel 
shock. The direction of the magnetic field determines the polarization (helicity) of 
the modes. The model parameters are kiVgo — 100, the particle distribution index 
a = 4.0, and the plasma parameter [3 = ag/v'^ = 1.0. The solid and dot-dashed 
curves show the two unstable modes. The normalizing parameter 70 is determined 
by Eq. (|13|l . See Bykov et al. (2011) for a comparison with the Bohm diffusion case 
r? = l 



I2OIOL ICastro and Sland[2010l) . This emission is likely du e to some combina- 
tion of pion-decay an d relativistic brcmsstrahlung (e.g.. ISykov et aLll2000l 
lUchivama eralll2010[) . 

Charge exchange reactions in the dovi^nstream region, from ne utrals cross- 
i ng th e shock and interacting, could result in a fast, cold ion beam. lOhira et all 
( 2009f ) presented a linear analysis of collisionless plasma instabilities between 
the cold beam and the hot downstream plasma. They found that, under SNR 
conditions, either the resonant instability or the Weibel-type instability are 
growing. The authors concluded that the mechanism may amplify the down- 
stream magnetic field to more than 100 /xG, changing the shock structure 
and the synchrotron spectra and profiles. However, as we have pointed out 
the turbulence amplification in the shock precursor is the main factor deter- 
mining the maximum CR energy and the efficiency of energy conversion to 
accelerated particles. 



3 Magnetic field amplification in diffusive shock acceleration 

The plasma instabilities in particle acceleration sources discussed above may 
result in magnetic turbulence and strong magnetic field amplification. The 
high efficiency of the acceleration means that nonlinear models are needed 
to predict the statistical characteristics of the amplified magnetic fields and 
their observational appearances. Diffusive shock acceleration models in su- 
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Fig. 2 The figure, from Bykov et al. (2010), shows growtli rates of the paral- 
lel propagating modes as a function of the wavenumber to illustrate the effect of 
short-scale turbulence on the long-wavelength instabihty. The model parameters are 
fco^'gO = 100 and a = 4.0. The solid and dot-dashed curves are simulated for two 
modes in the model with short-scale turbulence of ^ = 5 and = 10 to demonstrate 
the behavior of the modes in the intermediate regime. For comparison, the dashed 
and dotted curves are calculated fo r the mod el without the short-scale turbulence, 
i.e., with ATfl = and -> CX3 (c.f.. lBeiil[200l) . 



pernova remnants imply wide dynamic ranges for the fluctuations and the 
particle momentum spectra extending for more than four decades. These wide 
dynamic ranges present a severe problem for the direct use of particle-in cell 
plasma simulations and some simplifications arc required. 

The nonlinear modeling of the short-scale Bell instability was done using^ 



MHD simulations (e.g.. iBcll,2004 . Zirak ashvili and Ptuskini,2008., Zirakashvili et al.l 

assuming a fixed CR current as an external parameter. The MHD 
approach is justified since the CR current has only a weak response to 
the short-scale fluctuations. These models studied the spectral evolution 
in the short-scale range as well as the transformation of the turbulence 
through the subshock. The evolution of the Bell modes downstream from 
the shock was addressed b y Pelletier et al.l ( 2006() . iMarcowith et al.l ( 2006() . 



and lMarcowith and Cassi (|2010| ) 

A nonlinear Monte Carlo model of strong forward shocks in young su- 
pernova remnants with efficient particle acceleration, where a nonresonant 
instability driven by the CR curre nt amplifies magnetic tu rbulence in the 
shock precursor, was developed bv IVladimirov et al] (|2009( ). Particle injec- 
tion, magnetic field amplification, and the nonlinear feedback of particles 
and fields on the bul k flow were derived consistent ly in the Monte Carlo 
model as presented by IVladimirov et al] (|2006l l2008l ) . It was found that the 
shock structure depends critically on the efficiency of turbulence cascading. 
If cascading is suppressed, MFA is strong, the shock precursor is stratified, 
and the turbulence spectrum contains several discrete peaks. In Fig. [31 sim- 
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Fig. 3 Turbulence spectra, kW, at different locations relative to the subshock 
position as simulated by Vladimirov et al. (2009). The far upstream seed turbulence 
(at all k) is at the level indicated by the horizontal lines. 



ulated turbulence spectra, kW, are shown at different locations relative to 
the subshock position. 

In the solid curve, cascading is fully suppressed. In the dotted curve, 
the cascading from large to small scales is efficient and has the differential 
form corresponding to Kolmogorov's model. In the case with no cascading, 
the dissipation was assumed to be zero. With cascading, viscous dissipation 
was assumed. The peaks, apparent in the model with suppressed cascading, 
should influence synchrotron X-ray images of SNR shells, allowing observa- 
tional tests of cascading and other assumptions intrinsic to the nonlinear 
model of nonresonant wave growth. In the next section, we discuss some 
models of X-ray synchrotron imaging of SNR shells. 

Tw o opposite limits of cascading were considered in the model bv I Vladimirov et al.l 
(|2009f ) because our knowledge of turbulent cascading in coUisionless plasmas 
is still very limited. It is limited, in part, because of the rather narrow dy- 
namical range of the available nonlinear simulations. Information from direct 
measurements in the interplanetary medium can help constrain models of 
turbulent cascading, but this is also limited. The cascading of plasma wave 
spectra upstream of int erplanetary shocks d riven by coronal mass ejection 
events were analyzed bv lBamert et al.l (200l). They studied the competition 
between the upstream wave generation by suprathermal protons accelerated 
at the shock, and the cascading of wa ve energy in the inertial range of solar 
wind turbulence. iBamert et al.l (|2008[ ) concluded that amplified solar wind 
turbulence upstream of interplanetary traveling shocks is better described by 
a Iroshnikov-Kraichnan-type cascade rather than the Kolmogorov-type. 
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The turbulent cascade is exp ected to be anisotropic, as is the case in the 
interplanetary medium (see e.g.. lMatthaeus et alJll99fil l Horburv et alJl2008L 
|Podestall2009l . and the references therein). Turbulent cascade dissipation will 
heat ions and electrons both in the upstream and downstream regions of 
the shock. Since the magnetic field amplification region in the shock precur- 
sor is of a finite size ~ {c/vs)A, the upstream plasma temperature depends 
strongly on the heating rate and the plasma adyection time. The cascade 
rate depends on the field amplitude, growing mode polarization, and 
the plasma parameters. If t~^{c/v1)A < 1, the upstream plasma heating is 
inefficient. However, in the opposite limit with strong upstream turbulence, 
the plasma ion heating via an appropriate Landau resonance can be strong 
if the plasma parameter /? = ag/i;^ > I. Electron he ating is dominant in the 
highly mag netized upstream plasma if /3 < 1 fee.. fOuataert and Gruzinovl 
I1999L IHowcs ,2008,) . The M onte-Carlo simulat i ons o f the nonlinear structure 
of CR modified shocks by IVladimirov et al] (|2008l) demonstrated that the 
CR acceleration efficiency and the shock modification are sensitive to ion 
heating by the upstream turbulence. X-ray line emission from shocked shell s 
of SNRs can be used to constrain the DSA models fe.g.. lEllison et aDl201Ci l. 

There is an emerging class of SNRs, such as SN1006, RXJ1713. 72-3946, 
Vela Jr, and others, that are dominated by non-thermal emission in X-rays, 
probably of synchrotron origin. Such emission results from electrons and/or 
positrons that arc accelerated to well above TeV energies. Models of super- 
nova shells show that these leptons radiate in strong magnetic fields and the 
radiation losses produce a marked cut off in the lepton spectrum. Figure [3] 
(from .Ellison e t al, _20 10[) shows the X-ray emission from RXJ1713.72-3946 
as observed by the Suzaku spacecraft. This figure, and the models given in 
lEUison et al.l ( 2010t ). indicate that the lack of thermal X-ray emission lines 
can be important for determining the origin of the non-thermal emission and 
the strength of the magnetic fields in the SNR. Since the amplified magnetic 
fields in DSA must be fluctuating, we next discuss the effect of fluctuations 
on the synchrotron images of SNRs. 



4 Synchrotron images, spectra and polarization in supernova 
shells 

Non-thermal emission in many sources of interest such as SNRs, gamma-ray 
bursts, and AGNs originates from synchrotron radiation of ultra-relativistic 
particles in turbulent magnetic fields. The effect of a random magnetic field 
on synchrotron images, emission spectra, and polarization can be very sub- 
stantial. 



4.1 Synchrotron mapping of twinkling supernova shells 

It is instructive to illustrate first the effect of intermittcncy caused by the 
random magnetic field in terms of a power-law electron spectrum with spec- 
tral index F. The synchrotron emissivity /(r, t, v) ct shows that 
the local emissivity is relatively very high for large B and large F (a large F 
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Fig. 4 The top two panels show fits to the Suzaku SNR J1713 observations with 
a hadronic model for both Coulomb and instant temperature equilibration but 
ignoring the X-line emission. The blue (heavy wt.) solid curve is the synchrotron 
continuum, the black solid curve is the thermal emission (only lines above 10"'* MeV 
are included), and the dotted curve is the underlying bremsstrahlung continuum. 
The observed emission would be the sum (not shown in the top two panels) of 
the solid black and blue curves. The bottom panel shows the leptonic model (with 
Coulomb equilibration) where parameters have been chosen to be consistent with 
the Suzaku observations. For the hadronic model, the radiation intensity is multi- 
plied by 0.95 to match the observations. For the leptonic model, a normalization 
factor of 0.2 is required to match the observations. The Suzaku data have been 
adjusted for interstellar extin ction so no extinction is applied to the model in this 
plot. See lEUison eFall (pOTg i for full details. 



implies a steep spectrum). Therefore, the high-order statistical moments of 
the field dominate the synchrotron emissivity in the spectral cut-off regime 
and it is possible for a single strong local field maximum to produce a feature 
(dot or clump) that stands out on the map even after integrating the local 
emissivity over the line of sight. In lower energy maps, the contribution of 
a single maximum can be smoothed or washed out by contributions from 
a number of weaker field maxima integrated over the line of sight. A high- 
energy map, however, can be highly intermittent because the synchrotron 
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Fig. 5 Emission in the vicinity of a bright synchrotron dot Dl at different photon 
energies (see Bykov et al. 2008 for details). The four panels show, from left to right, 
■ J(Rx, v) maps of synchrotron emission at 0.5, 5, 20 and 50 keV, respectively. 



emissivity depends on high-order moments of the random magnetic field at 
the cut-off frequency regime. The electron/positron spectra in the cut-off 
regime are usually exponential, and the exact shapes of the spectra domi- 
nated by synchrotron/Compton losses depend on the diffusion coefficient. A 
power-law approximation can only apply for a narrow electron energy range 
in the cut-off regime, where the effective spectral index F is increasing with 

the electron energy. 

In o rder to calculate synchrotron emission maps and spectra. iBvkov et al.l 
modeled a system of finite size, filled with a random magnetic field 



and electro n spectra simulated with a DSA model accounting for synchrotron/Compton 
losses (see iBvkov et al.l |2000[) . A forward shock of spherical geometry was 
assumed in the image construction. The random field was composed of a su- 
perposition of magnetic fluctuations (plane waves) with random phases and a 
given spectrum of amplitudes. Non-steady, localized structures (dots, clumps 
and filaments), in which the magnetic field reaches exceptionally high val- 
ues, typically arise in the statistically stationary random field sample. These 
non-steady magnetic field concentrations dominate the synchrotron emission 
(integrated along the line of sight) from the highest energy electrons in the 
cut-off regime of the electron distribution, resulting in an intermittent, twin- 
kling, clumpy appearance even for a steady distribution of ultra-relativistic 
electrons. The spectral and temporal evolution of these high intensity syn- 
chrotron events, and the synchrotron images they produce, vary drastically 
at different photon energies. 

In Figure [5] the synchrotron emission of a small region around a dot 
labeled Dl is shown for emission at 0.5, 5, 20 and 50 keV. There are ob- 
servable differences in the synchrotron maps indicating that some features 
are bright at high energies and much less prominent at lower energies and 
vice versa. The physical reason for this difference comes from the fact that in 
the high-energy cut-off region of the electron spectrum, the local synchrotron 
emissivity depends strongly on the local field value as where the 

effective index F is increasing with the photon energy in the cut-off spectral 
regime. 

The simulated light curves of X-ray clumps may explain tho s e obse rved 
in X-ray images of some su pernova remnants bv lUchivama et al.l ( 20071 ) and 



iPatnaude and FesenI (j2008h . It is important to note that if the clumps are 
due to magnetic field intcrmittcncy, the time scale of their variation is not 
reflecting the ultra-relativistic electron loss time. The distinct characteris- 
tic of the modeled synchrotron emission is that its strong intermittency re- 
sults directly from the exceptionally high magnetic field amplifications that 



19 



randomly occur. The peaks in synchrotron emission maps occur due to high- 
order moments of the magnetic field probability distribution function (PDF). 
Even if the PDF of projections of the local magnetic field are nearly Gaus- 
sian, the corresponding PDF of synchrotron peaks, simulated for a spatially 
homogeneous relativistic particle distribution, has strong departures from 
the Gaussian at large intensity amplitudes. This is because of the nonlinear 
dependence of the synchrotron emissivity on the local magnetic field in the 
high-energy cut-off regime of the electron spectrum. Intermittency of this 
nature, with very different appearances, is rathe r a common phenomenon in 

stochastic media (see e.g., IZel'dovich etinilQSTl). 

The models bv iBvkov erall (|2008b[) and lBvkov et al.l (|2009| ) make a few 
basic predictions. One is that the overall efficiency of synchrotron radiation 
from the cut-off regime in the electron spectrum can be strongly enhanced 
in a turbulent field with some value of ^J{B^, compared to emission from a 

uniform field, Bq, where So = \/ {B^)- The second is that strong variations in 
the brightness of small structures can occur on time scales much shorter than 
variations in the underlying particle distribution. The variability time scale is 
shorter for higher energy synchrotron images. An estimate of the time scales 
of these intensity variations is consistent with the rapid time variability seen 
in some young SNRs by Chandra. The strong energy dependence predicted 
may be important for the future missions NuSTAR, Astro-H and others that 
will image SNR shells up to 50 keV. 



4.2 Synchrotron X-ray Polarization in DSA 

Synchrotron radiation in a regular, quasi-hom ogeneous magnetic field is po- 
larized (e.g., iGinzburg and Svrovatskiilll965[ ). It has long been known that 
random directions of magnetic fields, in addition to Faraday rotation, may 
stron gly reduce the average polarization of synchrotron emission source s 
re.g.. lWestfoldll959llCrusius and Schhckeiseill986llStroman and Pohlll2009D . 



This explains the relatively low polarization frequ ently observed for radio 
synchrotron sources. However, as it was shown by iBvkov et al ] (I2009D . the 



turbulent magnetic fields that reduce the average polarization can result 
in highly polarized, intermittent, patchy structures potentially observable in 
high resolution X-ray images. In terms of a power-law electron spectrum with 
spectral index F, the degree of polarization is given by 77 w {r + l)/{r+7/3) 
and, therefore, U increases with F. Thus the synchrotron images will be 
highly polarized at high energies. 

In order to construct maps of polarized synchrotron emission from SNR 
shel ls, it is convenient to use the loca l densities of the Stokes parameters (see 
e.g.. lGinzburg and Svrovatskiilll965l ). Because of the additive property of the 
Stokes parameters /, Q, U , V for incoherent photons, we can integrate these 
weighted with the distribution function of radiating particles over the line of 
sight across the source. The degree of polarization is then determined in a 
standard way as U = \/ -I- t/^ -I- II. Starting from the simulated ran- 
dom magnetic-field. iBvkov et all ( 20091 ) have constructed maps of polarized 
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Fig. 6 Simulated maps of polarized synchrotron emission in a random magnetic 
field at 5 keV. Intensity, ■ 1(11^,1, u), is shown with a linear color scale in the left 
panel. The central panel shows the product of intensity and polarization degree. 
The right panel shows the degree of polarization indicated by the color bar. The 
stochastic magnetic field sample has ^J{IP) = 3 x 10~^ G and spectral index of 
magnetic fluctuations is 5 = f.O (from Bykov et al. (2009)). 



X-ray emission of SNR shells. These are highly elumpy with polarizations up 
to 50% for energetic >TeV electrons in the cut-off regime. 

The detection of polarization in X-ray sources as faint as 1 milliCrab is an 
aim of the Gravity and Extreme Magnetism Small Explorer ( GEMS) mission, 
recently selected for flight in 2014 by NASA. It will make a sensitive search 
for X-ray polarization using thin foil mirrors and Time Projection Chambe r 
polarimeters with high efficiency in the 2-10 keV band ( Kallman et al.|[2Q10t ). 
The polarimeters under consideration, like XPOL aboard the planned IXO 
missio i{3 and POLARIX proposed to measure X-ray polarization with an 
angular resolution of about 20 arcsec in a field of view of 15 x 15 arcmin 
with the minimum detectable p olarization of 12% for a source of 1 mCrab 
within 100 ks of observing time ( Costa et al.ll2010[) . have good perspectives 
to detect polarized synchrotron emission from SNR shells. 

The intermittent appearance of the polarized X-ray emission maps of 
young SNR shells can be studied in detail with imagers of a few arcsecond 
resolution, though even arcmin rcsolutioi i images can provid e important in- 
formation as is illustrated in Fig. [5] (see iBvkov et al.ll2009() . The polarized 
emission clumps of arcsecond scales are time variable on a year or longer 
(depending on the observed photon energy, magnetic field amplification fac- 
tor, and the plasma density in the shell) allowing for rather long exposures 
even in the hard X-ray energy band. Hard X-ray observations in the spectral 
cut-off regime are the most informative to study the magnetic fluctuation 
spectra and the acceleration mechanisms of ultra-relativistic particles. 



^ http;//ixo. gsfc.nasa.gov/science/performanceRequirements. html 
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5 Magnetic field amplification in star forming galaxies 

The origin and evolution of maKnetic field structures in galaxies is an impor- 

tant u nsolved problem fe.g.. lRuzmaikin et al ■ll988llReesll2006LlKulsrud and Zweibel 



|2008[ ). Radio observations of the diffuse polar ized radio e mission from the 
disks of some spiral galaxies, as reviewed by iBeckl ( 2008( ). revealed large- 



scale coherent magnetic patterns that could be generated by dynamos, and 
dynamo models have been developed to explain the main observational fea- 
tures of the global magnetic fields of spiral galaxies (e.g., Ruzmaikin et al.l 
ll988llBeck et al.|[l996llBrandenburg and Subramanianll2005 ). In most galax- 



ies, however, the field structure is more complicated. For example, radio 
synchrotron observations of nearby galaxies reveal dynamically important 
magnetic fields of 10-30 /.iG total strength in the spiral arms. These spiral 
arm fields, however, tend to have random orientations, while ordered fields 
(observed in radio polarization) are strongest in interarm regions and follow 
the orientation of the adjacent gas spiral arms. 

A variety of observations (see [Wi drow 2 002i , for a review) suggest that 
magnetic fields are present in all galaxies including ellipticals and in galaxy 
clusters. Furthermore, it remains to be understood, in the context of the 
dynamo models, why well organized fields of surprisingly high strengths are 
observed in normal galaxies when the Universe was much younger than its 
present age (see e.g.. lBernet et al.ll2008llWolfe et al.ll2008( ). 



5.1 Magnetic fields in irregular galaxies 

Another issue to be addressed concerns the magnetic fields of dwarf irregular 
gal axies. Studies of magnetic fields in nearby dwarfs, as discussed recently 
bv IChvzvl ( 201Cl[) . revealed typically weak fields with the mean value of the 



total field strength about three times smaller than in normal spirals. The 
slow rotation rate in many irregular galaxies makes standard dynamo-type 
amplification of the magnetic field unlikely. On the other hand, dwarfs with 
stronger fields are associated with violent star-forming activity and these tend 
to be more massive and evolved systems. Their magnetic fields are thought 
to be reg ulated mainly by the surface density of the galactic star-formation 
rate (e.g.. IChvzvl[20Tflh . 

Star forming activity wit h clu stered supernova explosions in superbubbles 
observed in gala xies (see e . g. . | Heilesii 1990i) may affect the mean magnetic field 



as suggested by Ferriere (Il992 ). More over, numerical magnetohydrodynam- 
ical models bv ISieikowski et ahl ( 2010f ) of cosmic-ray driven dynamos in the 
interstellar medium of an irregular galaxy indicated the possibility of mag- 
netic field amplification. The process of cosmic-ray acceleration may play 
a role in generating galactic magnetic fields at different stages of cosmic 
evolution. Non-resonant mechanisms of magnetic field amplification allow 
field enhancements on scales much larger than the gyroradii of the ultra- 
relativistic particles accelerated in the sources. Such fields amplified by the 
CR current and CR pressure anisotropy will have scales well above a parsec 
in the star forming regions, and well above a kiloparsec in clusters of galaxies. 
Observations of supernova remnants and galactic CRs suggest that there is 
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a high efBciency of conversion of the kinetic power of supernova shocks into 
CRs. Therefore, the specific mechanisms discussed above may convert a siz- 
able fraction of the power released by SNRs into intermediate scale magnetic 
fields that are much larger than stellar scale fields. Some large-scale dynamo 
mechanism may also be necessary to produce the largest scale coherent fields 
seen in spiral galaxies. 



5.2 Magnetic fields in the first generation Pop III supernova remnants 



Simulations of n onlinear cosmic structure formation bv lAbel et al.l ( 20021) and 



IXu et al.l (I2008D showed that the high-redshift analog of a molecular cloud 



is able to produce individual massive primordial stars, offering a natural ex- 
planation for the absence of purely metal free low mass stars in the Milky 
Way. The production of small-scale dynamos by turbulence created by accre- 
tion shocks during gravitational collapse may amplify initially weak magnetic 
fields in the p rotostellar cloud in th e formation process of the first stars and 
galaxies fi.e.. ISchleicher et al.ll2010[ ). Massive stars exploding as supernovae 
provide energy, momentum, entropy, and metals at the very earliest stages 
of galactic evolution. Observations of high-redshift gamma-ray bursts likely 
indicate that relativistic particles and magnetic fields are produced during 
the collapse of massive primordial stars. 

An important problem is whether there is a way to transfer a significant 
fraction of the high kinetic energy of the supernova ejecta into the magnetic 
fields of the remnant. The questions to be addressed are the following. 

— (i) Can collisionless shocks be created by the supersonic ejecta in the very 
weakly magnetized plasma of /3 >> 1? 

— (ii) Can relativistic particles be accelerated in Pop III supernova rem- 
nants? 

— (iii) Are there effective mechanisms of magnetic field amplification in Pop 
III supernova remnants? 

Tw o-dimensional particle-in-cell simulations performed bv lKato and Takabd 
(|2010[ ) to study weakly magnetized perpendicular shocks demonstrated the 



development of a collisionless shock with /3 ss 26 (and an Alfven Mach num- 
ber of about 100). They showed that current filaments form in the foot region 
of the shock due to the ion-beam- Weibel instability (or the ion filamenta- 
tion instability) and these filaments generated a strong magnetic field there. 
Strong fluctuating field generation is necessary in the standard paradigm of 
a collisionless shock since the only way to decelerate and thermalize the cold 
upstream flow is through the electromagnetic fields. In the downstream re- 
gion, these current filaments also generated a tangled magnetic field that was 
about 15 ti mes stronger than the u pstream magnetic field. The shock sim- 
ulations bv iKato and Takabd ( 20101 ) reproduced well the Rankine-Hugoniot 



relations and indicated a substantial damping of the strong Weibel-type field 
fluctuations in the shock ramp (in the case of /? ~ 26, the magnetic field con- 
tribution in the Rankine-Hugoniot relations is small). The authors noticed 
that a fraction of the ions are accelerated slightly on reflection at the shock. 
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If the temporal and spatial scales of the simulation were increased, it is 
expected that this fraction of reflected superthermal particles would continue 
to be accelerated by the Fermi-type DSA process. In this case, the magnetic 
field would continue to be amplified in a symbiotic relationship where the 
magnetic turbulence required to scatter and acceler ate the CRs is created 
by the accelerated CRs themselves (e.s.. 16^1,2004. ,Vladimirov et ahlBoOSl 
iBvkov et al.ll201l| ). If this amplification operates as expected, a sizable frac- 
tion of the energy released by the supernova explosion can be transferred to 
energetic particles even with the initially weak magnetic field fluctuations of 
Bq ^ 0.1-1 nG expected in the protostellar cloud. 

The largest scale magnetic fluctuations amplifled by the long wavelength 
non-resonant instabilities (that are discussed above in SJH) can approach the 
scale of the shock. Therefore, the supernova remnant and its vicinity, with 
a scale size of about 30 pc, can be flUed with a magnetic field of strength 
~ 0.1 /iG generated by instabilities driven by CR acceleration. The longest 
scale fields generated by the non-resonant CR instabilities in the shock pre- 
cursor are in the coUisional plasma regime since their wavelengths are longer 
than the Coulomb mean free path of the ambient thermal plasma. The long- 
wavelength fields should have, by many orders of magnitude, longer damping 
times than the short-lived Weibel-type fluctuations that are located just in 
a very narrow vicinity of the shock jump. The magnetic flelds amplified in 
the SNR shock precursor can be considered as mesoscale fields since their 
scales are much larger than the stellar field sizes and larger than the CR gy- 
roradii, but are below the observed largest scale coherent galactic magnetic 
fields in spiral galaxies presumably generated by galactic dynamos or by a 
global galactic CR current. The mesoscale fields can be the dominant fields 
in irregular galaxies and can be spread by the galactic winds through galactic 
groups and clusters. Furthermore, CR-driven winds in high-rcdshift galaxies 
may have important implications for the metal and magnetic field enrichment 
of the intergalactic medium (see e.g.. lSamui et al ] |2010[lDubois and TevBsi^l 
120101) . In contrast, magnetic fields much lower than expected from CR insta- 
bilities were predicted in MHD simulations of the Bierma nn battery mech- 
anism in the first generation of supernova remnants (i.e., iHanavama et ahl 
I2009D . 

Magnetic field amplification by CRs produced by the first supernova ex- 
plosions during the re -ionization epoch of the universe was discussed by 
(|Miniati and Bellll2010[) . The CRs escape into the intergalactic medium car- 



rying an electric current that results in magnetic field amplification. The 
authors found that magnetic fields are robustly generated throughout inter- 
galactic space at a rate of 10"^'' — 10~^^ Gauss/Gyr, until the temperature of 
the intergalacti c medium is raise d by cosmic re-ionization. In another shock- 
related process. IRvu et all ( 20081 ) discussed a scenario in which turbulent-flow 



motions are induced via the cascade of vorticity generated at cosmological 
shocks during the formation of large-scale structure. It was argued that this 
turbulent dynamo model might provide average magnetic fleld strengths of 
a few microgauss inside clusters and groups, approximately 0.1 fiG around 
clusters and groups, and approximately 10 nG in filaments. 
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6 Discussion 

Magnetic fields, both regular and stochastic, play a key role in charged parti- 
cle accelerators providing both particle confinement and energy gain. The 
enormous energy release of gravitational (AGN, GRB, SNRs, clusters of 
galaxies, etc.) or thermonuclear origin (Type la SNRs, etc.), produces col- 
lisionless shocks and accelerates particles in many sources. It has long been 
known that strong shocks can only form in coUisionlcss plasmas if magnetic 
turbulence is generated and some thermal particles are injected and acceler- 
ated to superthermal energies. The strong interrelation between shock for- 
mation and structure, magnetic turbulence, and energetic particles makes 
understanding the self-consistent production of magnetic turbulence in dif- 
fusive shock acceleration both essential and difficult. 

While a great deal of progress has been made in this field, important 
challenges remain. In this brief review, we have discussed some recent work 
including resonant streaming and non-rcsonant CR-currcnt driven instabil- 
ities. An important example of this is the Bell mechanism that efficiently 
amplifies short-scale turbulence by the CR current. Magnetic turbulence am- 
plification may govern the nonlinear structure of CR-dominated flows, heat 
the upstream plasma, enhance synchrotron radiation, and limit high-energy 
electron energies. Furthermore, this mechanism might operate hand-in-hand 
with an instability for producing long-scale turbulence. The long-scale turbu- 
lence can be responsible for the highest energy particles accelerated by DSA. 
The basic DSA prediction that the highest energy accelerated ions can have 
a hard spectrum implies that the highest energy CRs may contain a size- 
able fraction of the shock ram pressure. This implies that a consistent model 
of magnetic fields in the highly supersonic flow with accelerated particles is 
nonlinear and multi-scale. 

The fleld of nonlinear shock acceleration is evolving rapidly with several 
prominent areas of active study. The long-wavelength instabilities discussed 
in Section [2] are yet to be combined with the resonant and short- wavelength 
instabilities in a nonlinear shock simulation. Such a combination will deter- 
mine the relative importance of each and how the overall shock structure is 
influenced by them. The instabilities produced by escaping cosmic rays are 
only just now starting to be explored. Only with a consistent model combin- 
ing the turbulence produced by escaping cosmic rays on the escaping cosmic 
rays themselves, can a reliable estimate of the maximum particle energy a 
given shock can produce be obtained. 

In principle, the plasma physics of shock structure, magnetic turbulence, 
particle injection and acceleration, indeed, all aspects collisionless plasmas, 
can be determined with particle-in-cell (PIC) simulations. A great deal has 
been done in this field and we have mentioned a small fraction of this work 
in Section [^31 The computational difficulties in producing PIC simulations 
with parameters close enough to systems such as SNRs are severe and we see 
future progress coming from synthesis of PIC models with semi-analytic and 
Monte Carlo techniques. 

We have also discussed some recent work on the observational conse- 
quences of magnetic turbulence and synchrotron emission. The important 



25 



result is that fluctuating magnetic fields strongly modify the spectra of syn- 
chrotron emission in the case of steep electron spectra that are typical for 
X-ray emitting, TeV particles in young SNRs. X-ray images of synchrotron 
dominated SNR shells are predicted to be highly intermittent with polarized 
twinkling structures of different scales if the strong amplification of magnetic 
turbulence is indeed operating in DSA. 

The current generation of X-ray telescopes - Chandra with its superb an- 
gular resolution and XMM-Newton with its high sensitivity ~ has uncovered 
a new world of synchrotron structures in young SNR shells. Further progress 
is certain to come from improved polarization and timing observations of 
SNRs at X-ray energies where the Faraday rotation is not essential and char- 
acteristic time scales are short enough. Observations at GeV-TeV energies 
are becoming more sensitive and theories of the highest energy cosmic rays 
must be improved accordingly. In all aspects of theory and observations, the 
turbulent magnetic field is the complicated glue tying things together. 
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